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Abstract 

In this paper we tackle the ANOVA problem for directional data (with particular emphasis 
on geological data) by having recourse to the Le Cam methodology usually reserved for linear 
multivariate analysis. We construct locally and asymptotically most stringent parametric tests 
for ANOVA for directional data within the class of rotationally symmetric distributions. We turn 
these parametric tests into semi-parametric ones by (i) using a studentization argument (which 
leads to what we call pseudo-FvML tests) and by (ii) resorting to the invariance principle (which 
leads to efflcient rank-based tests). Within each construction the semi-parametric tests inherit 
optimality under a given distribution (the FvML distribution in the first case, any rotationally 
symmetric distribution in the second) from their parametric antecedents and also improve on 
the latter by being valid under the whole class of rotationally symmetric distributions. Asymp- 
totic relative efficiencies are calculated and the finite-sample behavior of the proposed tests is 
investigated by means of a Monte Carlo simulation. We conclude by applying our findings on a 
real-data example involving geological data. 

Keywords : Directional statistics, local asymptotic normality, pseudo-FvML tests, rank-based 
inference, ANOVA. 

1 Introduction 

Spherical or directional data naturally arise in a broad range of earth sciences such as geology (see, 
e.g., Watson 1983 or Fisher and Hall 1990), astrophysics, meteorology, oceanography or studies of 
animal behavior (see, e.g., Merrificld 2006 and the references provided therein) or even in neuro- 
science (see Leong and Carlile 1998). Although primitive statistical analysis of directional data can 
already be traced back to early 19th century works by the likes of C. F. Gauss and D. Bernoulli, 
the methodical and systematic study of such non-linear data by means of tools tailored for their 
specificities only begun in the 1950s under the impetus of Sir Ronald Fisher's pioneering work (see 
Fisher 1953). We refer the reader to the monographs Fisher et al. (1987) and Mardia and Jupp (2000) 
for a thorough introduction and comprehensive overview of this discipline. 

An important area of application of spherical statistics is in geology (for instance for the study 
of palaeomagnetic data, see McFadden and Jones 1981 or the more recent Acton 2011) wherein the 
data are usually modeled as realizations of random vectors X taking values on the surface of the 
unit hypersphere S''^^ := {v G M'^ : v'v = 1}, the distribution of X depending only on its angular 
distance from a fixed point 9 G S''^^ which is to be viewed as a "north pole" for the problem under 
study. A natural, flexible and realistic family of probability distributions for such data is the class 
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of so-called rotationally symmetric distributions introduced by Saw (1978) ~ see Section [2] below for 
definitions and notations. Roughly speaking such distributions allow to model all spherical data that 
are spread out uniformly around a central parameter with the concentration of the data waning 
as the angular distance from the north pole increases. Within this setup, an important question 
goes as follows : "do several measurements of remanent magnetization come from a same source of 
magnetism?" More precisely, suppose that there are m different data sets spread around i sources 
of magnetism 0i £ S^~^ , i = 1, . . . ,m. The question then becomes that of testing for the problem 
Ho . 01 = . . . = m against "Hi : 3 \ < i ^ j < m such that 0i ^ 0j, that is, an AN OVA problem for 
directional data. 

This important problem has, obviously, already been considered in the literature (see Mardia 
and Jupp 2000, chapter 10, for an overview). The difficulty of the task, however, entails that most 
available methods are either of parametric nature or suffer from computational difhculties/slowness 
such as WcUncr (1979) 's permutation test or Beran and Fisher (1998) 's bootstrap test. To the best 
of our knowledge, the only computationally simple and asymptotically distribution-free test for the 
general null hypothesis Ho above is the test given in Watson (1983). The purpose of the present 
paper is to complement this literature by constructing tests that are optimal under a given m-tuple 
of distributions-(Pi , . . . ,Pm) say-but remain valid (in the sense that they meet the nominal level 
constraint) under the general null hypothesis Ho involving a large family of spherical distributions. 
In particular, the tests we propose are asymptotically distribution-free within the semi-parametric 
class of rotationally symmetric distributions. Obviously the applicability of our ANOVA procedures 
is not reserved to geological data alone, but directly extends to any type of directional data for which 
the assumption of rotational symmetry with location parameter seems to be reasonable. 

The backbone of our approach is the so-called Le Cam methodology (see Le Cam 1986), as 
adapted to the spherical setup by Ley et al. (2013). Of utmost importance for our aims here is the 
uniform local asymptotic normality (ULAN) of a sequence of rotationally symmetric distributions 
established therein and which we adapt to our present purpose in Section |3l In the same Section [3] we 
also adapt results from Hallin et al. (2010) to determine the general form of a so-called asymptotically 
most stringent parametric test for the above hypothesis scheme Ho against Hi- Due to its parametric 
nature the optimality of the (Pi, . . . , P,„) -parametric test is thwarted by its non- validity under any 
m-tuple (Qi, . . . , Qm) distinct from (Pi, . . . , P,„). In order to palliate this problem we have recourse 
to two classical tools which we adapt to the spherical setting: first a studentization argument, which 
leads to so-called pseudo-Fisher-von Mis es- Lang evin (pseudo-FvML) tests, and second the invariance 
principle, yielding optimal rank-based tests. Both families of tests are of semi-parametric nature. 

The idea behind the pseudo-FvML test has the same flavor as the pseudo- Gaussian tests in 
the classical "linear" framework (see, for instance, Muirhead and Waternaux 1980 or Hallin and 
Paindaveine 2008 for more information on pseudo-Gaussian procedures) . More concretely, since the 
FvML distribution is generally considered as the spherical analogue of the Gaussian distribution (see 
Section [2] for an explanation), our first approach consists in using the FvML as basis distribution 
and "correcting" the (parametric) FvML most stringent test, optimal under a m-tuple (Pi, . . . ,Pm) 
of FvML distributions, in such a way that the resulting test 0^"^ remains valid under the entire class 
of rotationally symmetric distributions. We obtain the asymptotic distribution of the asymptotically 
most stringent pseudo-FvML test statistic Q'"^ under the null and under contiguous alternatives. 
As it turns out, the test statistic Q^"^ and the test statistic Qwitson provided in Watson (1983) are 
asymptotically equivalent under the null (and therefore under contiguous alternatives). As a direct 
consequence, we hereby obtain, in passing, that Watson (1983)'s test (/>witson ^^^^ enjoys the property 
of being asymptotically most stringent in the FvML case. 

The optimality property of 0watson ^^"^ '/'^"^ ^y construction, restricted to situations in which 
the underlying m-tuple of distributions is FvML. In the sequel we make use of the well-known 
invariance principle to construct a more flexible family of test statistics. To this end we first obtain 
a group of monotone transformations which generates the null hypothesis. Then we construct tests 
based on the maximal invariant associated with this group. The resulting tests (that are based on 
spherical signs and ranks) are, similarly as (j^wltson ^^'^ asymptotically valid under any m- 
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tuple of rotationally symmetric densities. Our approach here, however, further entails that for any 
given m-tuple (Pi, . . . , Pm) of rotationally symmetric distributions (not necessarily FvML ones) it 
suffices to choose the appropriate m-tuple = (Ki, . . . ,Km) of score functions to guarantee that 
the resulting test is asymptotically most stringent under (Pi, . . . , Pm)- 

The rest of the paper is organized as follows. In Section [51 we define the class of rotationally 
symmetric distributions and collect the main assumptions of the paper. In Section [3l we summarize 
asymptotic results in the context of rotationally symmetric distributions and show how to construct 
the announced optimal parametric tests for the ANOVA problem. We then extend the latter to 
pseudo-FvML tests in Section |3] and to rank-based tests in Section [5J and study their respective 
asymptotic behavior in each section. Asymptotic relative efficiencies are provided in Section [S) The 
theoretical results are corroborated via a Monte Carlo simulation in Section[6l A real data application 
is considered in Section [71 Finally an appendix collects the proofs. 



2 Rotational symmetry 

Throughout, the m(> 2) samples of data points X,;i, . . . , Xi„;, i = 1, . . . , m, are assumed to belong 
to the unit sphere S''^^ of R*^, k > 2, and to satisfy 

Assumption A. (Rotational symmetry) For alH = 1, . . . , m, X^, • . • , Xi„; are i.i.d. with com- 
mon distribution Pg.;/. characterized by a density (with respect to the usual surface area measure 
on spheres) 

^^ckj, m^'e,), xe^'^-i, (2.1) 

where 0i G S'^^^ is a location parameter and /i : [— 1, 1] — > is absolutely continuous and (strictly) 
monotone increasing. Then, if X has density ()2.ip . the density of X'^; is of the form 

where ujk = /T{k/2) is the surface area of iS'^"^ and B{-,-) is the beta function. The corre- 

sponding cumulative distribution function (cdf) is denoted by Fi{t) , i = 1, m. 

The functions fi are called angular functions (because the distribution of each X^ depends 
only on the angle between it and the location 6i G S''~^). Throughout the rest of this paper, we 
denote by J"™ the collection of m-tuples of angular functions / := (/i, /2, . ■ • , fm)- Although not 
necessary for the definition to make sense, monotonicity of fi ensures that surface areas in the vicinity 
of the location parameter 6i are allocated a higher probability mass than more remote regions of 
the sphere. This property happens to be very appealing from the modeling point of view. The 
assumption of rotational symmetry also entails appealing stochastic properties. Indeed, as shown in 
Watson (1983), for a random vector X distributed according to some Psi ji as in Assumption A, not 
only is the multivariate sign vector Sei(X) := (X — (X'0i)0i)/||X — {'X.'6i)6i\ \ uniformly distributed 
on iS*i" := {v G K*^ I || v|| = 1, v'6i = 0} but also the angular distance X'^i and the sign vector Sg. (X) 
are stochastically independent. 

The class of rotationally symmetric distributions contains a wide variety of useful spherical dis- 
tributions including the wrapped normal distribution, the FvML, the linear, the logarithmic and the 
logistic (a definition of the latter three is provided in Section [5] below) . The most popular and most 
used rotationally symmetric distribution is the aforementioned FvML distribution (named, according 
to Watson 1983, after von Mises 1918, Fisher 1953, and Langevin 1905), whose density is of the form 

/fvML(k)(x;0) = Ck{K)cxp{Kx.'0), X e S''"'^, 

where k > is a concentration or dispersion parameter, £ S''~^ a location parameter and Ck{n) 
is the corresponding normalizing constant. For ease of reference we shall, in what follows, rather 
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use the notation 1/)^ instead of /fvML(k)- This choice of notation is motivated both by the wish for 
notational simphcity but also serves to further underhne the analogy between the FvML distribution 
as a spherical and the Gaussian distribution as a linear distribution. This analogy is mainly due 
to the fact that the FvML distribution is the only spherical distribution for which the spherical 
empirical mean ^Moan := X)r=i -^i/ll Sr=i (based on observations Xi,...,X„ G 5'^"^) is the 
Maximum Likelihood Estimator (MLE) of its spherical location parameter, similarly as the Gaussian 
distribution is the only (linear) distribution in which the empirical mean Sr=i (based on ob- 
servations Xi, . . . , X„ G M'") is the MLE for the (linear) location parameter. We refer the interested 
reader to Breitenberger (1963), Bingham and Mardia (1975) or Duerinckx and Ley (2013) for details 
and references on this topic; see also Schaeben (1992) for a discussion on spherical analogues of the 
Gaussian distribution. 



3 ULAN and optimal parametric tests 

Throughout this paper a test (jj* is called optimal if it is most stringent for testing Hq against Hi 
within the class of tests Ca of level a, that is if 

sup r^.(P) < sup r0(P) V(/. G C^, (3.2) 
PG-Hi Pe-Hi 

where r^^ (P) stands for the regret of the test (j)o under P G "Hi defined as r^,, (P) := [sup^^,^^ Ep [(/)]] — 
Ep[(/)o], the deficiency in power of (po under P compared to the highest possible (for tests belonging 
to Ca) power under P. 

As stated in the Introduction, the main ingredient for the construction of optimal (in the sense 
of p.2p ) parametric tests for the null hypothesis Hq . Oi = 62 = ■ ■ ■ = Om consists in establishing 
the ULAN property of the parametric model 

({p^|,J^ig5'=-i},...,{p£,J0™g5^-^}) 

for a fixed m-tuple of (possibly different) angular functions / = (/i,...,/m), where P^"/; stands 
for the joint distribution of X^i, . . . ,Xi„;, i = 1, . . . ,to, for a fixed m-tuple of angular functions 
(/i, . . . , fm)- Letting •& := {9[, . . . ,^'„)', we further denote by P^"j the joint law combining Pg^^.f^, 

(n) . ~ 

. . . , Pg .J .In order to be able to state our results, we need to impose a certain amount of control 
on the respective sample sizes n^, i = 1, . . . , m. This we achieve via the following 

Assumption B. Letting n = jy^i '^i' f^^' all i = 1, . . . , m the ratio r^-"-* := n^/n converges to a finite 
constant as n — )■ cw. 



In particular Assumption B entails that the specific sizes Ui arc, up to a point, irrelevant; hence 
in what precedes and in what follows, wc simply use the superscript '■"^ for the different quantities 
at play and do not specify whether they are associated with a given n^. In the sequel we let 
diag(Ai, . . . , A,„) stand for the m x m block-diagonal matrix with blocks Ai, . . . , A,„, and use the 
notation diag((r("))-i/2l,, . . . , (rtY'/^h)- 

Informally, a sequence of rotationally symmetric models |P^"j I S (5'^~^)'"| is ULAN if, uni- 
formly in t?^") = (0^"^', . . .,0l^^'y G (5'=-!)™ such that 1?^"^ - ■& = 0(71-1/2), ^j^^ log-likelihood 

allows a specific form of (probabilistic) Taylor expansion (see equation p. 41) below) as a function 
of t(") := (tj"^', . . . jtm^'y G R™'^. Of course the local perturbations t^") must be chosen so that 



4 



i?'"-* + n i/2j/(")t(") remains on (5*^ and thus, in particular, the tl""* need to satisfy 

= (^1"' + n-^/^t("))'(0m n-^/^ti" V 1 
^2r^;'^'{0^^Yt!r^+n-'(t^yi^ (3.3) 

for all i ~ l,...,m. Consequently, t^"'' must be such that 2n,^^^^(0,["-')'t^"'' + o{n~^^'^) = : for 
^i"^ + i/^tj-"'' to remain in S'^^^, the perturbation t|"^ must belong, up to a o{n^ ^^^) quantity, to 
the tangent space to S'^^^ at 0^"'. 

The domain of the parameter being the non- linear manifold (5*^^^) it is all but easy to establish 
the ULAN property of a sequence of rotationally symmetric models. A natural way to handle this 
difficulty consists, as in Ley et al. (2013), in resorting to a re-parameterization of the problem in 
terms of spherical coordinates r/. say, for which it is possible to prove ULAN, subject to the following 
technical condition on the angular functions. 

Assumption C. The Fisher information associated with the spherical location parameter is finite; 
this finiteness is ensured if, for i = 1, . . . ,m and letting Lpj^ := jij ji (fi is the a.e. -derivative of /,), 

After obtaining the ULAN property for the 77-parameterization, one can use a lemma from Hallin et 
al. (2010) to transpose the ULAN property in the spherical r/-coordinates back in terms of the original 
^-coordinates. Finally the inner-sample independence and the mutual independence between the m 
samples entail that we can deduce the required ULAN property which is relevant for our purposes 
(this we state without proof because it follows directly from Proposition 2.2 of Ley et al. 2013). 

Proposition 3.1 Let Assumptions A, B and C hold. Then the model |P^"j | G (5'^'~^)™| is 
ULAN with central sequence A^"| :~ ^(A^"^^^)', . . . , {Ag^\^ )'^ , where 

:-'^7'/'E^/.(X^.^')(l-(X^.^')')'^'S,.(X,,), i = l,...,m, 

and Fisher information matrix T^:/ := di&giVg-^. j-^^ . . . ^^e„^■.f,n) where 

More precisely, for any 

^(n) g (5fe-i)m g^^f^ that ■d^'^^ ^ ■& = 0(7i-i/2) and any bounded sequences 
t(") = (tj"^', . . .,tin^'y as m {23), we have 

108 ( = (t<...).A«,, - l(tW)'r,,tl". + „p(l), (3.4) 

where A^"!^,^ A- J\fmkiO,T^.j), both under P^"|, as n ^ 00. 

Proposition 13.11 provides us with all the necessary tools for building optimal /-parametric pro- 
cedures (i.e. under any m-tuple of densities with respective specified angular functions fi, . . . , /,„) 
for testing Hq . Oi = ... = 9m against Hi : 31 < i ^ j < to such that 7^ Oj. Intuitively, 
this follows from the fact that the second-order expansion of the log-likelihood ratio for the model 
|p^"j I G (iS'^^^)™! strongly resembles the log-likelihood ratio for the classical Gaussian shift ex- 
periment, for which optimal procedures are well-known and are based on the corresponding first-order 
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term. Now clearly the null hypothesis Hq is the intersection between {S^ and the linear subspace 
(of M'"*^) 

C := {v = (v;,...,v;„)' I vi,...,v„ e R^-'and vj = . . . = v,„} Ifc) 

where we put 1,„ := (1, ■ • ■ , 1)' G M™, 7W(A) for the linear subspace spanned by the columns of 
the matrix A and A (X) B for the Kronecker product between A and B. Such a restriction, namely 
an intersection between a linear subspace and a non-linear manifold, has already been considered in 
Hallin et al. (2010) in the context of Principal Component Analysis (in that paper, the authors ob- 
tained very general results related to hypothesis testing in ULAN families with curved experiments). 
In particular from their results we can deduce that, in order to obtain a locally and asymptotically 
most stringent test in the present context, one has to consider the locally and asymptotically most 
stringent test for the (linear) null hypothesis defined by the intersection between C and the tangent 
to (S''~^)™'. Let denote the common value of 0i, . . . ,6m under the null. In the vicinity of !,„ (E) 0, 
the intersection between C and the tangent to (iS*^^^)"* is given by 

[{0' + ,i-i/2(^("))-i/2t(")', ,,,^e' + ?.-i/2(rL"V/'t(;:)')', (3.5) 

e't("^ = . . . = e'ti?) = 0, (r("V'/'t("' = . . . = (ri?V'/'t,l:^}. 

Solving the system p.5p yields 

^(«)t(«) = ((r^"V'/'t^"^...,(rH)-l/2t(;^)')' eA4(i,„® (i,.-ee'))- (3.6) 

Loosely speaking we have "transcripted" the initial null hypothesis Hq into a linear restriction of the 
form p.6p in terms of local perturbations t'^"\ for which Le Cam's asymptotic theory then provides a 
locally and asymptotically optimal parametric test under fixed /. Using Proposition 13.11 and letting 

J-m ®{Ik-00') and T^"^ := (i/("))-iT tf, an asymptotically most stringent test 0/ is then 
obtained by rejecting Ho as soon as (A~ stands for the Moore-Penrose pseudo-inverse of A) 

A;^, (t^.^^ - Tg ((Tg)'r,^,Tg) " (Xg)') A.^; (3.7) 

exceeds the a-upper quantile of a chi-square distribution with (to — l)(fc — 1) degrees of freedom. 
Hence the optimal parametric tests are now known. 

There nevertheless remains much work to do. Indeed not only does the optimality of our test 
4>f only hold under the TO-tuple of angular densities / = (/i, . . . , /m), but also this parametric test 
suffers from the (severe) drawback of being only valid under that pre-specified m-tuple. Since it 
is highly unrealistic in practice to assume that the underlying densities are known, these tests are 
useless for practitioners. Moreover, we so far have assumed known the common value of the spherical 
location under the null, which is unrealistic, too. The next two sections contain two distinct solutions 
allowing to set these problems right. 



4 Pseudo-FvML tests 

For a given TO-tuple of FvML densities ((/)„j , . 
Ki, . . . , Km > (where we do not assume ki ~ 



. . , ) with respective concentration parameters 
. . . — Km), the score functions (/s^^ reduce to the 
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constants k^, i ~ 1, . . . , m, and hence the central sequences for each sample take the simplified form 



rii 

=: K,i{lk-9iO[)n\^'^'Xi 

= Kt{Ik-Oi0[)ny^{'Xi~ei), i = l,...,m. 

Optimal FvML-bascd procedures (in the sense of p.2|) ) for Ho are then built upon A^"j_ ~ (Ag"^^ , • ■ • , ^'g^^'-^ )'j 
where ((/)ki 

Before proceeding we here again draw the reader's attention to the fact that a parametric test built 
upon A^"^ will only be valid under the m-tuple and becomes non- valid even if only the concentration 
parameters change. In this section, this non-validity problem will be overcome in the following way. 
We will first study the asymptotic behavior of A^"^ under any given m-tuple g = {gi, . . . ,gm) € J^"^ 

and consider the newly obtained quadratic form in A^"^. Clearly, this quadratic form will now 

depend on the asymptotic variance of A^"^ under g, hence again, for each g, we are confronted to 

an only- for- valid test statistic. The next step then consists in applying a studentization argument, 
meaning that we replace the asymptotic variance quantity by an appropriate estimator. We then 
study the asymptotic behavior of the new quadratic form under any m-tuple of rotationally symmetric 
distributions. As we will show, the final outcome of this procedure will be tests which happen to 
be optimal under any m-tuple of FvML distributions (that is, for any values ki, . . . , k„i > 0) and 
valid under the entire class of rotationally symmetric distributions; these tests are our so-called 
pseudo-FvML tests. 

For the sake of readability, we adopt in the sequel the notations E/[-] for expectation under 
the angular function / and ■do Im 'E) (where we recall that represents the common value 
of 6i, ... ,0m under the null). The following result characterizes, for a given m-tuple of angular 
functions g e J^'", the asymptotic properties of the FvML-based central sequence A^"^^, both under 

P^"^^ and P^"+„-i/2y(„)(.(,i).g with t^") as in p.3p for each sample. 

Proposition 4.1 Let Assumptions A, B and C hold. Then, letting Bk.g^ := 1 — Eg. [(X^^^)^] for 
i = 1, . . . ,m, we have that A^"''^ is 

(i) asymptotically normal under P^^^j.^ with mean zero and covariance matrix 

rs„;g :=diag {rig^,...r;.gj , 

B 

where T*g.g^ := (Ifc - 00'), i = 1, . . . ,m; 

(ii) asymptotically normal under P^"^_|_„-i/2y(n)(.(n).g (t^"^ as in \3.3\) ) with mean r^Q;0,gt (t := 

.(«) 



(t'l, . . . ,t',j)' with ti lim„_^oo , i ^ 1, . . . ,m) and covariance matrixT^^.g, where, putting 
^e)2)^,^(X^6 



Ck,g, := E,J(1 - {Xl^er)ipg^ (X'^M forz^l,...,m. 
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withTg.^^^^g. := '^'fc,"^^"' (Ifc -09'), i = l,...,m. 

Sec the Appendix for the proof. As the null hypothesis only specifics that the spherical locations 
coincide, we need to estimate the unknown common value 6. Therefore, we assume the existence of 
an estimator oi such that the following assumption holds. 

Assumption D. The estimator ■& = 1,„ (g) ^, with e S''~^, is n'^/'^(i/^"'^)~'^ -consistent: for all 
■00 = Im ®0 ena, - t9o) = Op(l), as n ^ oo under P^"^^ for any g e J"™. 

Typical examples of estimators satisfying Assumption D belong to the class of Af -estimators (see 
Chang 2004) or i?-estimators (see Ley et al. 2013). Put simply, instead of A^"-}^ we have to work with 

A^"^ for some estimator 'd satisfying Assumption D. The next crucial result quantifies in how far 

this replacement affects the asymptotic properties established in Proposition l4.1l fa proof is provided 
in the Appendix). 

Proposition 4.2 Let Assumptions A, B and C hold and let'd — !,„ ®0 he an estimator of'0Q such 
that Assumption D holds. Then 

(i) letting T'-"'' :— r^"^ : . . . : \J Vm^ I^^ , ■^tf"-^ satisfies, under P.^"!g and as n ^ oo, 

Ai^^ -Ai"\ = -r^ T("'V^f0-0Kop(l), 

where 

rl.,:=diag(r^-,...,rfc) 

with r^.;; K,Eg^ [X;./] (I^ - eO'), i^l,...,m; 

It) for all d e {S'^-^r, T^,^^^ = = T*^,^. 

Following the inspiration of Hallin and Paindaveine (2008) (where a very general theory for 
pseudo- Gaussian procedures is described) we are in a position to use Proposition 14.21 to construct 
our pseudo-FvML tests. To this end define, for i = 1, . . . , m, the quantities Ek^g^ := Eg. [X-^-fl] , and 

set, for notational simplicity, D^.g^ :~ Ek.gjBk.g^ and H^^g := X^I^i ^i"^^fc,gi^fc,Si ■ Then, letting 
— fr* \-r— t(") [i't^") Vf— fr* \-r— t^"^ ]-fyi"-) vf— fr* '\- 

\'-^o;g) '-'&o-g'-So;vi\'-aoM' '■■Oo;g^ '^o;g' ■0o;g ■6o;vi y^^o;^' ■So;g^^o;gJ ' 
the g-valid test statistic for "Ho : 0i = . . . = 0m we propose is the quadratic form 

Q(")(g) := (A^"^)'ri^ aI"^ 

It is easy to verify that Q*-"' {g) does not depend explicitly on the underlying concentrations ki , . . . , 
but still depends on the quantities Bk^g^ and Ek,g. , i ~ I, . . . ,m. This obviously hampers the validity 
of the statistic outside of g. The last step thus consists in estimating these quantities. Consistent (via 

the Law of Large Numbers) estimators for each of them are provided by Bk^g^ := l — n~^ 'J2^Li0^ijO)^ 

and Ek^gi n^^ IZ^jLii^ijO) i i ~ 1, . . . , m. For the sake of readability, we naturally also use the 

notations bk^g^ := Ek,gjBk,g,, i = l,...,m, and i?0,g := J2iLif'f^^blg_Bk,g,. Putting X^ := 
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rij ^ X]j=i for all 2 = 1, . . . , m, straightforward calculations then show that our pscudo-FvML 
test statistic for the m-sample spherical location problem is 

m A rn A A 

which no more depends on g. 

The following proposition, whose proof is given in the Appendix, finally yields the asymptotic 
properties of this quadratic form under the entire class of rotationally symmetric distributions, 
showing that the test is well valid under that broad set of distributions. 

Proposition 4.3 Let Assumptions A, B and C hold and let i? be an estimator of jJq such that 
Assumption D holds. Then 

(i) Q^"' is asymptotically chi-square with (m — l)(fc — 1) degrees of freedom under 

(ii) Q^"^ is asymptotically non-central chi-square with fc — 1 degrees of freedom and non-centrality 
parameter 

under P^^^J^^_-^^2i,{n)^{n).g; where t^^'^ is as in hS.S^) and t := lim„_i.oo t^"^/ 

(Hi) the test 0^"' which rejects the null hypothesis as soon as Q^"-* exceeds the a-upper quantile of 
the chi-square distribution with (m — l)(fc— 1) degrees of freedom has asymptotic level a under 

(iv) 0*^"' is locally and asymptotically most stringent, at asymptotic level a, 
Utf„e«o Uge.F-'iP^oLj' '^ga^^st alternatives of the form [j^^^„{P'^'^^}- 

Remark 1 It is easy to verify that Q'"^ is asymptotically equivalent (the difference is a op(l) quan- 
tity) to the test statistic for the same problem proposed in Watson (1983) under the null (and there- 
fore also under contiguous alternatives). Thus, although the construction we propose is different, our 
pseudo-FvML tests coincide with Watson's proposal. In passing, we have therefore also proved the 
asymptotic most stringency of the latter. 



5 Rank- based tests 

The pseudo-FvML test constructed in the previous section is valid under any m-tuplc of (non- 
necessarily equal) rotationally symmetric distributions and retains the optimality properties of the 
FvML most stringent parametric test in the FvML case. Although the FvML assumption is often 
reasonable in practice, our aim in the present section is to depart from this assumption and provide 
tests which are optimal under any distribution. 

We start from any given m-tuplc / G J^'" and our objective is to turn the /-parametric tests 
into tests which are still valid under any m-tuple of (non-necessarily equal) rotationally symmetric 
distributions and which remain optimal under /. To obtain such a test, we have recourse here to 
the second of the aforementioned tools to turn our parametric tests into semi-parametric ones: the 
invariance principle. This principle advocates that, if the sub-model identified by the null hypothesis 
is invariant under the action of a group of transformations Qt, one should exclusively use procedures 
whose outcome does not change along the orbits of that group Qt- This is the case if and only if 
these procedures are measurable with respect to the maximal invariant associated with Qt- The 
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invariance principle is accompanied by an appealing corollary for our purposes here: provided that 
the group Gt is a generating group for Ho, the invariant procedures are distribution- free under the 
nuU. 

Invariance with respect to "common rotations" is crucial in this context. More precisely, letting 
O e SOk := {A € M'^^'^jA'A = 1^, det(A) = 1}, the nuU hypothesis is unquestionably invariant 
with respect to a transformation of the form 

However, this group is not a generating for Ho as it docs not take into account the underlying angular 
functions /. which are an infinite-dimensional nuisance under Ho- This group is actually rather 

generating for Ui^qG-Ho ^i"-/ ^^^^ fixed /. Now, denote as in the previous section the common value 

oi0i,...,e^ under the nuU as 0. Then X.^ = (X^0)6I + ^1 - (X^6')2Se(X,j) for all j = 1, ... , n, 

and i = 1, . . . , m. Let Qh ih. '■= {hi, . . . , ft-m)) be the group of transformations of the form 

5,. : X,, ^ gn^ (X,, ) = h, (X^0)0 + y^l - (/i,(X^0))2Se (X,, ) , * = 1 , . . . , m, 

where the hi : [—1, 1] — > [—1, 1] are monotone continuous nondecrcasing functions such that hi{l) = 1 
and hi{—l) = —1. For any 77i-tuple of (possibly different) transformations (gh^^, ■ ■ ■ , gh,„) G Gh, 
it is easy to verify that ||gft.(Xij)|| = 1; thus, gh^ is a monotone transformation from S'^^^ to 
S'^~^ , i = l,...,m. Note furthermore that gi^ does not modify the signs Se(Xy). Hence the 
group of transformations Gh is a generating group for IJ f^jrm P^"''/ and the null is invariant un- 
der the action of Gh- Letting Rij denote the rank of ^ijO among X'j^d, . . . , XJ„ .0, i = 1, . . . , m, 
it is now easy to conclude that the maximal invariant associated with Gh is the vector of signs 
Se(Xii), . . . , Se(Xi„j), . . . , Se(Xmi), ■ . ■ , Se(Xm„^) and ranks . . . , Rim, ■ ■ ■ , Rmi, ■ ■ ■ , Rmn^- 

As a consequence, we choose to base our tests in this section on a rank-based version of the 
central sequence A^"^^, namely on 

with 

4 a. = r^;'^" E (^) S,(X.,), ^ = 1, . . . , m, 

j=i \ ' / 

where := {Ki, . . . , Km) is a m-tuplc of score (generating) functions satisfying 

Assumption E. The score functions Ki, i = 1, . . . ,m, are continuous functions from [0, 1] to M. 

The following result, which is a direct corollary (using again the inner-sample independence and the 
mutual independence between the m samples) of Proposition 3.1 in Ley et al. (2013), characterizes 
the asymptotic behavior of A^"^^ under any m-tuple of densities with respective angular functions 

gii---,grn- 

Proposition 5.1 Let Assumptions A, B, C and E hold and consider g — {gi, . . . ,gm) G J-™ ■ Then 
the rank-based central sequence A^"^^ 

(i) is such that A^"^^^ — A^"-'^^.^ = op(l) under P.^"^^ as n ^ oo, where (Gi standing for the 
common cdf of the ^ijO 's under P^J.g, i = 1, • • • , fn) 
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with 

rii 

i=i 

In particular, for K ^ Kf := {Kf,, . . . , Kf^) with Kf^{u) := (pf^{Fr^ {u)){l - {Pr^ {u))^y/^ , 
A^"'^ is asymptotically equivalent to the efficient central sequence A^"^^ under P^J.j:. 

(a) is asymptotically normal under P^"^^ with mean zero and covariance matrix 



, m. 



T^,.,K ■■= diag {^^{ik 90% ee': 

where Jk{Ki) := Kf{u)du. 

(Hi) is asymptotically normal w?i(ier P^"^_j_^_j^^2^(„j^(„)^ (t^"^ as in i3.3\} } with meanT^^-x^gt (t 
lim„_i.oo t^") ) and covariance matrix 

r..;^, -.^ diag - 09'), - 99') 

where Jk{Ki, gi) := K.i{u)Kg. {u)du for i = 1, . . . , m. 
(iv) satisfies, under P'^J,^ as n ^ oo, the asymptotic linearity property 

4i:V„-V..(.,t(".;K - 4£k = -r.o;K,,t(") + 0P(1), 
/ort(") = (t("^',...,t(;^'')' as 11 



Similarly as for the pseudo-FvML test, our rank-based procedures are not complete since we 
still need to estimate the common value 6 of 9i, ... ,6m under Ho- To this end we will assume the 
existence of an estimator i? satisfying the following strengthened version of Assumption D : 

Assumption D'. Besides ri^/^(i/("^) ^-consistency under P^"^^ for any g g J"™, the estimator ■d S 
j-^fe-i-jm f^j-|;i^gj- locally and asymptotically discrete, meaning that it only takes a bounded number 
of distinct values in i?o-centcrcd balls of the form {t G M™'' : n^/^H (i/("')^^(t - i?o)ll < c}- 



Estimators satisfying the above assumption are easy to construct. Indeed the consistency is not 
a problem and the discretization condition is a purely technical requirement (needed to deal with 
these rank-based test statistics, see pages 125 and 188 of Lc Cam and Yang 2000 for a discussion) 
with little practical implications (in fixed-n practice, such discretizations arc irrelevant as the radius 
can be taken arbitrarily large). We will therefore tacitly assume that 6 G S'^~^ (and therefore 
d = l„i (X) 9) is locally and asymptotically discrete throughout this section. Following Lemma 4.4 
in Kreiss (1987), the local discreteness allows to replace in Part (iv) of Proposition 15.11 non-random 
perturbations of the form ■& + n^i/ V(")t(") with t(") such that ■d + n-i/2i/(")t(") still belongs to Ho 
by a n^/^(i/^"-')~^-consistent estimator ^ := 1,„ (g) 0. Based on the asymptotic result of Proposition 
15.11 and letting 
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the g-valid rank-based test statistic we propose for the present ANOVA problem corresponds to the 
quadratic form 

This test statistic stiU depends on the cross-information quantities 

Jk{Kf,,g^),...,MKf^,g,,) (5.8) 

and hence is only valid under fixed g. Therefore, exactly as for the pseudo-FvML tests of the previous 
section, the final step in our construction consists in estimating these quantities consistently. For 
this define, for any p > 0, 

Up) + nT^'^p [k - - ee') a * = l, . . . , m. 

Then, letting 9i{p) := 0i{p)/\\6i{p)\\, we consider the piecewise continuous quadratic form 

Consistent estimators of the quantities J';7"^(i^i,gi), . . . , J^^{Km, gm) (and therefore readily of (|5.8|) ') 
can be obtained by taking 

p, := mi{p > : h'f'\p) < 0} 

for i = l,...,m (see also Ley et al. 2013 for more details). Denoting by Jk{Ki,gi), for i = 

l,...,m, the resulting estimators, setting Hji.g := ^^'^ letting 

Ki (jiij/{ni + 1)^ Sg(Xy ), i — 1, ... ,171, {Rij naturally stands for the rank of X^^^ among X^j^^, . . . , X'j„.^), 

(n) 

the proposed rank test 0}^ rejects the null hypothesis of homogeneity of the locations when 
Qin) n n, ^.^-i ^ ri,njJk{Kug^)Ju{Kj,gj) 



exceeds the a-upper quantile of the chi-square distribution with (to — l)(fc — 1) degrees of freedom 



This asymptotic behavior under the null as well as the asymptotic distribution of Q under a 



sequence of contiguous alternatives are summarized in the following proposition. 

Proposition 5.2 Let Assumptions A, B, C and E hold and let'd be an estimator such that Assump- 
tion D' holds. Then 

(i) Q is asymptotically chi-square with (m — l)(fc — 1) degrees of freedom under 

(ii) Q^^^ is asymptotically non-central chi-square, still with (to — l)(fc — 1) degrees of freedom, hut 
with non-centrality parameter 

Mnrfer P^"^^^_i/2^(„)t(ii)-g; where t^"'^ is as in iS. 3\) andt := lim„_j.oo t^"-'; 

(Hi) the test which rejects the null hypothesis as soon as Q^'' exceeds the a-upper quantile of 
the chi-square distribution with (to — l)(fc — 1) degrees of freedom has asymptotic level a under 
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(iv) m particular, for Kj := (i^/i, . . ■ ,i^/,„) with Kf^{u) := (p f^{Fr\u)){l - {Ff^{u))'^)'^^'^, 

0^'' is locally and asymptotically most stringent, at asymptotic level a, for UiJog-Hq Ugej='™{-Pi"''g} 

';/J 



against alternatives of the form [J^^-u„{P^j}- 



Thanks to Proposition 15.11 the proof of this resuh fohows along the same Unes as that of Proposi- 
tion 14.31 and is therefore omitted. 

We conclude this section by comparing the optimal pseudo-FvML test 0'-"-' with optimal rank- 
based tests for several choices of / e J^™ by means of Pitman's asymptotic relative efficiency 

(ARE). Letting ARE^(,.g((^5."\ 4'2^'^) denote the ARE of a test 0^"^ with respect to another test "* 
under P^"-*, , we have that 

ARE^^.g(0^'^,0^"^) = l'do,t-K_f,g/ho,t:l,g- 

In the homogeneous case g = (gi, . . . , gi) (the angular density is the same for the m samples) and if 
the same score function — namely, Kf-^ — is used for the m rankings (the test is therefore denoted by 



,(") 



I, the ratio in ^ simplifies into 



ARE,„,(^(,"; /^H) ^ ^.^j-f^^-"^^ . (5.9) 

Numerical values of the AREs in (|5.9[) are reported in Table [1] for the three-dimensional setup under 
various angular densities and various choices of the score function Kf-^. More precisely, we consider 
the spherical linear, logarithmic and logistic distributions with respective angular functions 

/iin(a)(i) ■■=t + a, /log(a) (0 := log(t + a) and 

. . a exp(— & arccos(t)) 

/iogis(a,b)(i) := (i + aexp(-6arccos(t)))2- 

The constants a and b are chosen so that all the above functions are true angular functions satisfying 
Assumption A. The score functions associated with these angular functions are denoted by Kiin{a.) for 
/lin(a), ^iog(a) for /log(a) and Kiogis{a,h) for /iogis(a,b)- For the FvML distribution with concentration 
K, the score function will be denoted by K^^. 

Inspection of Table [T] confirms the theoretical results. As expected, the pseudo-FvML test i/)'^"' 
dominates the rank-based tests under FvML densities, whereas rank-based tests mostly outperform 
the pseudo-FvML test under other densities, especially so when they are based on the score function 
associated with the underlying density (in which case the rank-based tests are optimal). 



6 Simulation results 

In this section, we perform a Monte Carlo study to compare the small-sample behavior of the pseudo- 
FvML test and various rank-based tests 4> ^k' ^'-'^ '^^^ two-sample spherical location problem, that 
is, for an ANOVA with m = 2. For this purpose, we generated M — 2, 500 replications of four pairs of 
mutually independent samples (with respective sizes ni ~ 100 and n2 = 150) of (k =)3-dimcnsional 
rotationally symmetric random vectors 

ee-iji, ^=1,2,3,4, = 1, . . . , n^, i = l,2, 

with FvML densities and linear densities: the Si-ij-^'s have a FvML(15) distribution and the ei-2j2^s 
have a FvML(2) distribution; the e2;iji's have a Lin(2) distribution and the e2;2j2's have a Lin(l.l) 
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Tabic 1: Asymptotic relative efficiencies of (homogeneous) rank-based tests (j) ^'^ with respect to the 
pseudo-FvML test i/)^"^ under various three-dimensional rotationally symmetric densities. 



















Underlying density 








2 /f,i„(4) 




X ifloEi3(l,l) 


X -R'logiB(2,l) 


FvML(l) 


0.9744 


0.8787 


0.9813 


0.9979 


0.9027 


0.9321 


0.7364 


FvML(2) 


1 


0.9556 


0.9978 


0.9586 


0.9749 


0.9823 


0.8480 


FvML(6) 


0.9555 


1 


0.9381 


0.8517 


0.9768 


0.9911 


0.9280 


Lin(2) 


1.0539 


0.9909 


1.0562 


1.0215 


1.0212 


1.0247 


0.8796 


Lin(4) 


0.9709 


0.8627 


0.9795 


1.0128 


0.8856 


0.9231 


0.7097 


Log(2.5) 


1.1610 


1.1633 


1.1514 


1.0413 


1.1908 


1.1625 


1.0951 


Log(4) 


1.0182 


0.9216 


1.0261 


1.0347 


0.9503 


0.9741 


0.7851 


Logis(l,l) 


1.0768 


1.0865 


1.0635 


0.9991 


1.0701 


1.0962 


0.9778 


Logis(2,l) 


1.3182 


1.4426 


1.2946 


1.0893 


1.4294 


1.3865 


1.5544 



distribution; the e^-ij-^^s have a FvML(15) distribution and the £3:2^2 have a Lin(l.l) distribution 
and finally the e4:ijj's have a Lin(2) distribution and the £4:2^2 's have a FvML(2) distribution. 

The rotationally symmetric vectors et-^j- 's have all been generated with a common spherical 
location = (•\/3/2, 1/2, 0)'. Then, each replication of the Se-^ij.^s was transformed into 

f Xf;ijj = e^aji, ^ = 1,2,3,4, ji = l,...,ni 

I ^«;2j2;« = C)e£<?;2i2i ^ = 1,2,3,4, j2 = l,---,n2, ^ = 0,1,2,3, 

where 

/ cos«/16) -sin(7r^/16) 
= sin(7r$/16) cos«/16) 
\ 1 

Clearly, the spherical locations of the X£;ijj's and the Xf;2j2;o's coincide while the spherical loca- 
tion of the X£;2j2;c'S: C = li 2, 3, is different from the spherical location of the ^e-iji 's, characterizing 
alternatives to the null hypothesis of common spherical locations. Rejection frequencies based on 
the asymptotic chi-square critical values at nominal level 5% are reported in Table [5] below. The 
inspection of the latter reveals expected results: 

(i) The pseudo-FvML test and all the rank-based tests are valid under heterogeneous densities. 
They reach the 5% nominal level constraint under any considered pair of densities. 

(ii) The comparison of the empirical powers reveals that when based on scores associated with the 
underlying distributions, the rank-based test performs nicely. The pscudo-FvML test is clearly 
optimal in the FvML case. 



7 Real- data example 

In this section, we evaluate the usefulness of our tests on a real-data example. The data consist of 
measurements of remanent magnetization in red slits and claystones made at 2 different locations in 
Eastern New South Wales, Australia. These data have already been used in Embleton and McDonnell 
(1980). The rotationally symmetric assumption in the two samples seems to be appropriate since 
data are clearly concentrated. However, the specification of the angular functions is not reasonable. 

The main question for the practitioner is to test whether the remanent magnetization obtained 
in those samples comes from a single source of magnetism or not. Therefore, we test here the null 
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Table 2: Rejection frequencies (out oi M = 2,500 replications), under the null and under increas- 
ingly distant alternatives, of the pseudo-FvML test 0'"^ and various rank-based tests (f) ^ 

(based on FvML(15) and FvML(2) scores), ^l^^^^. ^^-j (based on Lin(2) and Lin(l.l) scores), 

4>^^ K \ (based on Lin(2) and FvML(2) scores) and </>J'^^ ^ (based on FvML(15) and 

Lin(l.l) scores). Sample sizes are ui = 100 and ?i2 = 150. 









Test 


True densities 


n 
u 


1 
1 


9 


Q 
O 


^(") 




.0592 


.2684 


.8052 


.9888 






.0696 


.2952 


.8276 


.9900 


^(i^Li„(2),-ffLi„(l.l)) 


( r cho 1 


.0536 


.2316 


.7660 


.9756 


/ (") 
2(ifL,„(2),iffe) 




.0656 


.2952 


.8160 


.9894 


1 (n) 

2(K^,,,Kl,„(1.1)) 




.0544 


.2308 


.7716 


.9772 






.0480 


.0596 


.0792 


.1312 






.0472 


.0568 


.0948 


.1340 


/ (n) 

2(KL.„(2),-ffL.„(l.l)) 


(Lin(2),Lin(l.l)) 


.0464 


.0604 


.0892 


.1424 


/ (") 

2(ifL.„(2),if,*2) 




.0520 


.0588 


.0920 


.1440 


1 (") 

2(i^*,r,.-'^L,„(l.l)) 




.0480 


.0580 


.0856 


.1340 


(^(") 




.0508 


.0684 


.1044 


.1512 


, (") 




.0540 


.0648 


.1012 


.1532 


1 (n) 

2(^fLi„(2),-ffLi„(l.l)) 


(Lin(2),</)2) 


.0512 


.0664 


.1084 


.1608 


1 (") 
2(^fLi„(2):if</,2) 




.0508 


.0656 


.1072 


.1620 


1 (") 

2(K*,,,,-ffL.„(l.l)) 




.0496 


.0628 


.1004 


.1516 






.0468 


.1008 


.2908 


.5760 


An) 




.0628 


.1288 


.3612 


.6788 


1 (n) 

2(i<:L.„(2),ifL.„(l.l)) 


(0i5,Lin(l.l)) 


.0512 


.1156 


.3636 


.6892 


1 (") 

2(ifL.„(2),iffe) 




.0616 


.1220 


.3620 


.6768 


1 (") 




.0504 


.1180 


.3660 


.6916 
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hypothesis Hq ■.61—62 agamst Hi -.61 ^62- For this purpose, we used the pseudo-FVML test t/)'^"^ 
and rank-based tests 1) fvML(io)) ^^'^ ^(nn(i 1) FvML(ioo)) based respectively on the couples 

of linear and FvML scores (lin(l.l), FvML(lO)) and (lin(l.l), FvML(lOO)). The corresponding test 
statistics are given by 

g(") = 5.96652, g Si(i.i),F,ML(io)) = 5-477525 and Q !ui(i.i),p.ML(ioo)) = 5-525854. 
At the asymptotic nominal level 5%, the tests 0^"-*, (iin(i.i),FvML(io)) !^ (Un(i.i),FvML(ioo)) 

do 

not reject the null hypothesis of equality of the modal directions since the 5%-upper quantile of the 
chi-square distribution with 2 degrees of freedom is equal to 5.991465. 



Magnetic remanence (data 1) Magnetic remanence (data 2) 




Figure 1: Measurements of remanent magnetization in red slits and claystones made at 2 different locations 
in Australia 

Appendix 

Proof of Proposition [47l] From Watson (1983) (and the beginning of Section [5]) we know that, 
under P^"'!^, the sign vectors Sg(X.ij) are independent of the scalar products X^-^^, Eg. [Se(Xij)] ~ 
and that 

E,, [S,(X,,)S^(X,,)] = ^(I,. - 66') 

for i = 1, . . . , m and for all j = 1, . . . , n^. These results readily allow to obtain Part (i) by applying 
the multivariate central limit theorem, while Part (ii) follows from the ULAN structure of the model 
in Proposition 13.11 and Le Cam's third Lemma. □ 

Proof of Proposition [472] We start by proving Part (i). First note that easy computations yield 
(for i = 1, . . . , m) 

i=i 

rii 
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where V^"' n^n-'J^^U [^^A «'^'(^ - ^nd wl"' 9 K,n-'(jy;u^^)r^/'{e 6). Now, 
combining the deha method (recah that 1^ — 06 is the Jacobian matrix of the mapping /i : M*^ — > 
X H- > evaluated at 0), the Law of Large Numbers and Shitsky's Lemma, we obtain that 



S 



k-l 



under Pi"L as n 



V 



(") 



K.E^JX:/] (Ife - 6e')nY\e -0)+ op(l) 



— >■ c«. Thus, the announced resuh follows as soon as we have shown that is 

op(l) under P^"'!^ as n ^ oo. Using the same arguments as for V^"'', we have under P^"^g and for 
n — cx) that 



6 I ^,nri ^(X^) (I, - 00') I n]/\6 - 0) + op(l) 



K^Eg, [y^l - (X',/)2(S,(X.y))'] n^^iO - 0) + op(l), 



which is op(l) from the boundedness of and since from Watson (1983) (see the proof of Proposi- 
tion 14.11 for more details) we know that 



E„ 



l-(X^0)2(Se(X,,))' 



E„ 



E,J(S,(X,,))']=0'. 



This concludes Part (i) of the proposition. Regarding Part (ii) , let X be a random vector distributed 
according to an FvML distribution with concentration k. Then, writing c for the normalization 
constant, a simple integration by parts yields 

Ck,^^ = kE^JI - {X'ef] = f^cj (l-M^) exp(Ku)(l-it2)('^-3)/2 



exp 



(Atu)(l-7.2)(fc-l)/2 du 



c(fc- 1) 



u exp 



(km)(1-u2)('=-3)/2 



(fc-l)E^jX'0]. 



The claim thus holds. 



□ 



Proof of Proposition 14.31 We start the proof by showing that the replacement of 6 with as 
well as the distinct estimators have no asymptotic cost on Q*-"-*. The consistency of Dk,gi, Ek,gi, 
i = 1, . . . , m, and H(f,.g together with the n^/^(i/'^"))~^-consistency of entail that, using Part (i) of 
Proposition 14.21 

g(") = {A^-l^^^Tl^J^-^V^^{0-O)y^ 
under P^"''^ as n — > oo. Now, standard algebra yields that 



17 



so that 

Q("H'?o)+op(l). 

Both resuhs from Proposition l4.1l entail that since r^^^.^r^^^.^ ^ is idempotent with trace (m— l)(fc— 1), 
(and therefore Q*^"^) is asymptotically chi-square with (m — l){k — 1) degrees of freedom 
under P^"''^, and asymptotically non-central chi-squarc, still with (m— l)(fc— 1) degrees of freedom, 

and with non-centrality parameter tT^^.^ gF^^ ^r^^.^ ^t under P^"^^_„-i/2t,(n)t(n) g- Parts (i) and (ii) 
follow. Now, Part (iii) is a direct consequence of Part (i). Part (ii) of Proposition 14.21 and simple 
computations yield that Q*^"' is asymptotically equivalent to the most stringent FvML test Q^"' in 
Part (iv) thus follows. □ 
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